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We investigate the presence of percolating states in disordered two-dimensional topological insula¬ 
tors. In particular, we uncover a close connection between these states and the so-called topological 
Anderson insulator (TAI), which is a topologically non-trivial phase induced by the presence of 
disorder. The decay of this phase could previously be connected to a delocalization of bulk states 
with increasing disorder strength. In this work we identify this delocalization to be the result of a 
percolation transition of states that circumnavigate the hills of the bulk disorder potential. 
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INTRODUCTION 

Topological insulators are recently discovered mate¬ 
rials with promising electronic properties'"^. The first 
physical system that was identified as a topological insu¬ 
lator was the HgTe/CdTe quantum well^ featuring one¬ 
dimensional edge states that are protected from backscat- 
tering by time-reversal symmetry, strongly stabilizing 
them against non-magnetic disorder. The existence of 
such conducting edge states is hallmarked by a quan¬ 
tized conductance plateau that has meanwhile been veri¬ 
fied experimentally.^ As discovered in a successive numer¬ 
ical study, ^ such edge states are not only immune from 
backscattering but can even be elicited by disorder in sys¬ 
tems that have no topologically distinct properties in the 
clean limit. This disorder-induced topological phase was 
first believed to be caused by Anderson localization, and 
was thus named ’’topological Anderson insulator” (TAI). 
A later study found, however, that the phase boundaries 
at the transition from an ordinary insulator to the TAI 
can be explained by an effective medium theory in which 
the presence of disorder leads to a re-normalization of the 
medium parameters.^ In this sense, the TAI appears due 
to a change of topology in the effective medium. 

While the transition from an ordinary insulator to the 
TAI could be explained by the aforementioned effective 
description,^ the transition from the TAI phase back to 
an ordinary insulating phase at very strong disorder val¬ 
ues proves more involved: the bulk states localize at 
intermediate disorder strength allowing for unimpeded 
edge-transport in the TAI phase, yet delocalize when 
disorder becomes even stronger.So far, the resulting 
breakdown of the TAI phase could be attributed to the 
coupling of counter-propagating edge states on oppos¬ 
ing edges through these delocalized bulk states, result¬ 
ing in a suppression of the edge states’ immunity from 
backscattering. This mechanism is responsible for an 
increased sensitivity to finite size effects^^d3 making the 
transition hard to explore numerically and leaving the 
true nature of this counter-intuitive delocalization un¬ 
clear. While first studies^d3,i4 interpreted the bulk de- 
localization as an intermediate metallic phase, a later 
study^^ considering larger systems pointed out that only 
a single bulk state is probably responsible for the delocal¬ 


ization and an intermediate metallic phase is not present. 
In addition, a spatially correlated potential and the as¬ 
sociated pronounced bulk delocalization turn out to de¬ 
stroy the TAI phase entirely.In this work, we resolve 
the puzzle associated with these different observations 
by identifying the emergence of percolating states as the 
origin of the delocalization and by clarifying the general 
connection between such states and the TAI phase. 

MODEL AND METHODS 

As a starting point we choose the well-studied dis¬ 
ordered HgTe/CdTe quantum well as described by the 
Bernevig-Hughes-Zhang (BHZ) model^ in terms of an ef¬ 
fective 4-band Hamiltonian 

= (I) 

with 

h{k) = le{k) + di{k)a^ (2) 

eik) =C-D{kl + kl) 

—*■ T 

di = ( Akx, -Aky, M{k) ) 

M{k) =m- B{kl + kl) 

and cr* representing the Pauli-matrices. Following 

Ref. 16, we choose the following set of realistic quan¬ 

tum well parameters in all our computations: A = 
364.5 nm • meV, B = —686.0 nm^ • meV and D = 
—512.0 nm^ • meV. The topology of the system is de¬ 
termined by the sign of the topological mass m: For 
m < 0 the bulk band gap of size 2|m| is topologically 
non-trivial and thus filled with gap-less edge states char¬ 
acterizing a two-dimensional topological insulator. On 
the other hand, for m > 0, the bulk band gap is topolog¬ 
ically trivial and does not contain any states leaving us 
with a system that is an ordinary insulator. 

Using the advanced modular recursive Green’s function 
method^we calculate the conductance G through 
two-dimensional rectangular ribbons of HgTe/CdTe 
quantum wells discretized on a square grid with dis¬ 
cretization constants Ax and Ay, width W = {nyPl)'Ay 
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and length L = • Ax. In accordance with previ¬ 

ous studies,the discretization constants are 
set to Ax = Ay = a = 5 nm. Two clean semi-infinite 
leads are attached to the left and right end of the ribbon. 
Following the Landauer-Biittiker formalism, the conduc¬ 
tance G in the limit of vanishing temperature is given 
by the total transmission T at the Fermi energy Ep. 
Our method also allows for a calculation of the scatter¬ 
ing wave functions 2 p{x,y) as well as the density of states 
p{E) = —Tr [Im G'^{x, x, E)] /tt where G^ is the retarded 
Green’s function. Disorder is modelled by static on-site 
energy values V{xi,yi) at each grid point {xi^yi) ran¬ 
domly chosen from the interval [—I//2,I//2] with U the 
disorder strength. In most studies,the values 
of V{xi^yi) are chosen without any spatial correlations 
between neighboring grid points [see Fig. la]. Here, we 
also consider spatial correlations in V{x,y) [see Fig. 1], 
characterized by a finite correlation length which can 
significantly affect the conduction properties. 


RESULTS 

In our simulations we first consider a quadratic region 
of width W = L = 1000 nm and an uncorrelated disorder 
potential (i.e., ^ 0). Two geometries will be stud¬ 

ied that only differ in their boundary conditions: a rib¬ 
bon for which hard wall boundary conditions along the 
edges are applied and a cylinder with periodic bound¬ 
ary conditions in ^-direction. A comparison between the 
disorder-averaged conductance (G) through these two ge¬ 
ometries has been used previously to distinguish between 
bulk and edge phenomena as the periodicity of the cylin¬ 
der eliminates the edges of the geometry. The results 
for (G) in a topological insulator with m = — lOmeV at 
Fermi energy Ep = 16meV are shown in Fig. 2(a) as a 
function of the disorder strength U for both geometries. 
The value of Ep is chosen such that for uncorrelated dis¬ 
order the TAI conductance plateau with (G) = le^/h 
clearly appears in the ribbon^’for disorder strength 
80meV < U < 280 meV [see red curve in Fig. 2(a)]. In 
the cylinder this plateau is clearly absent, since no edge 
states can exist in this edge-less geometry. While for dis¬ 
order values beyond this plateau the conductance drops 
monotonically in the ribbon, the conductance through 
the cylinder geometry [blue dashed curve in Fig. 2(a)] 
shows a renewed increase at the same disorder strength. 
This is the signature of the aforementioned bulk delocal¬ 
ization that has already been observed in uncorrelated 
potentials.A physical intuition for this transition is, 
however, still lacking, but will become clear by consid¬ 
ering disorder potentials that are spatially correlated [il¬ 
lustrations for uncorrelated and correlated disorder po¬ 
tentials are shown in Fig. 1]. 

In a previous work^^ we demonstrated that spatial cor¬ 
relation in the disorder potential can destroy the TAI 
phase entirely. Here we consider a situation, in which 
the correlations all but dissolve the plateau in the ribbon 


geometry [see the red curve in Fig. 2(b)]. For these pa¬ 
rameter values it is best visible that the dissolution of the 
plateau is accompanied by a delocalization of bulk states. 
As can be seen by comparing the blue dashed curves in 
Fig. 2(a) and (b), this delocalization happens at much 
lower disorder values for correlated potentials than for 
uncorrelated ones. In both cases, however, these delo¬ 
calized bulk states contribute to the conductance, but 
also suppress the edge conductance by coupling counter- 
propagating edge states with each other, thereby leading 
to a breakdown of the TAI conductance plateau. 

To get a better insight into this scenario, we now con¬ 
sider the scattering wave functions 2 p{xi,yi) during this 
delocalization transition [see Fig. 3(a) and (c) for two 
such states at 1/ = 220 meV]. A first optical inspection of 
these wave function images suggests that the associated 
flux is circumnavigating the hills of the underlying corre¬ 
lated potential,reminiscent of percolation states found 
in the Quantum Hall effectand in antidot topological 
insulator lattices.To make this first impression more 
quantitative, we analyze how the intensities of the wave 
functions 'ip{xi^yi) shown in Fig. 3 are correlated with the 
values V{xi,yi) of the underlying potential landscape. 
For this purpose we compute the weighted probability 
P{V) for a wave to encounter a potential value V with the 
weights of this probability distribution being given by the 
intensity \ 2 p{xi,yi)\‘^ of the wave function at a grid-point 
(xi^yi) with potential value V. The distribution P{V) 
resulting from an average over 1000 disorder realizations 
shows a surprisingly pronounced enhancement at posi¬ 
tive disorder values V approximately situated between 
Vmin ~ 25meV and Umax ~ 75meV [see Fig. 4], suggest¬ 
ing that disorder values from this interval give rise to 
clearly enhanced wave function intensities. Apparently 
the states responsible for the bulk delocalization tend to 
reside primarily at relatively high values of the disorder 
potential, i.e., in a certain altitude interval of the hills in 
the correlated potential landscape. Correspondingly, we 
find that the wave function intensities shown in Fig. 3(a) 
and (c) strongly resemble contour plots of the associated 
disorder potential, when we truncate that latter to the 
interval V G [25,75] meV [see Fig. 3(b) and (d)]. 

To identify the origin of this curious behavior, we first 
point to the fact that the above interval bounds, i.e., 
Umin ~ 25meV and Unax ~ 75meV, are astonishingly 
close to the minimal and maximal distances F^min = 
26meV and F^max = 78meV of the Fermi energy Ep = 
16meV to the energy range F^buik = [—10, —62] meV in 
which the valence band states are situated in the clean 
system [see the band structure of Fig. 5(a)]. This obser¬ 
vation suggests that the flux in our correlated potential is 
carried mostly by the disorder-analogues of these valence 
band states. Further evidence for this correspondence 
can be deduced when considering the rescaled probabil¬ 
ity distribution P{Ef — V), which measures, as above, 
the probability for a scattering state to reside at a poten¬ 
tial value U, but now shifted by the Fermi energy Ep. 
We find that this distribution, quite remarkably, stays 
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almost invariant with respect to a change of the Fermi 
energy [see a comparison between two different values of 
Ep in Fig. 5(c)]. This observation reflects the fact that a 
change oiEp just shifts the corresponding wave functions 
to different disorder values, but that the origin of states 
in the valence bands stays unchanged. Furthermore, the 
density of states p{E)^ shown in Fig. 5(b), and the distri¬ 
bution P{Ef — V), shown in Fig. 5(c), are very similar 
- even the small kinks in p{E) are clearly reproduced in 
P{Ef — V). Such kinks in the density of states are noth¬ 
ing else but van-Hove-singularities resulting from the flat 
bands in the band structure. We may thus conclude that 
the flat valence band states especially at the Brillouin 
zone (BZ) edges (where the maximum of P{Ef — V) is 
located) represent the most significant contribution to 
the intensity of the scattering wave functions. In addi¬ 
tion, we find that the distribution P{Ef — V) is not at 
all sensitive to the boundary conditions since it is almost 
exactly equal for the ribbon and for the cylinder (see 
Fig. 5(b)). The distribution P{Ef — V) thus turns out 
to be quite fundamental in that it has its origin in ba¬ 
sic system properties, which are given here by the band 
structure in Fig. 5(a) and by the flat band states con¬ 
tained in it. 

These observations allow us to construct a compre¬ 
hensive picture of the physics in the cylindrical geom¬ 
etry with a correlated potential [see blue dashed curve in 
Fig. 2(b)]: While in the clean limit pure bulk conduction 
is observed, the conductance drops down to a minimum 
at disorder strength oiU 60 meV due to the increasing 
localization of the bulk states [see Fig. 2(b)]. When at 
U ~ 100 meV the hills of the correlated potential are high 
enough to locally shift the Fermi energy into the valence 
band, the bulk is filled with localized states deriving from 
clean valence band bulk states spatially located around 
the hills of the underlying disorder potential. With grow¬ 
ing disorder strength, these localized states connect with 
each other and go through a percolation threshold, which 
is responsible for the delocalization transition and the 
increased bulk conductance. Only at very strong disor¬ 
der U ^ 300 meV the connection between these percolat¬ 
ing states weakens and the conductance again decreases. 
This picture is also strongly supported by previous stud¬ 
ies of the TAI in the uncorrelated case (see Refs. 23 and 
20): Considering the arithmetic and geometric average 
of the local density of states it was shown there that the 
states carrying the flux in the TAI are not single extended 
states throughout the whole TAI phase (as would be ex¬ 
pected for edge states) but for very strong disorder rather 
formed by clusters of well localized states. Our percolat¬ 
ing wave functions deriving from the valence band are 
perfect candidates for such linked, localized states. This 
picture is also corroborated by the flatness of the valence 
band states, which leads to the very small group velocity 
responsible for the wave function enhancements around 
the potential hills as seen in the examples of Fig. 3(a) 
and (c). 

The flatness of the states in the effective band structure 


is, in fact, also important for the theory put forward in 
the aforementioned studies:^^’^^ Considering disordered 
super-cell structures it was argued that flat and localized 
bands develop small sub-gaps that can be of topological 
non-trivial type. Hence these gaps have to be filled with 
edge states in the same way as the inverted band gap of a 
clean topological insulator is.^’^’^^ In this picture the TAI 
phase is thus characterized by edge states that appear in 
the energy gaps between localized bulk states and are 
consequently again immune from backscattering. 

At this point the question arises how the above results 
can be reconciled with our own model, which so far does 
not contain any reference to edge states in the percola¬ 
tion transition. To investigate this issue in detail, we per¬ 
formed additional calculations for a system where no edge 
effects can be present due to a topological mass, which 
we choose to take on the positive value of m = 10.0 meV. 
As shown in Fig. 2(c) this sign change of m significantly 
modifies the conductance properties. While previously 
for m < 0 and moderate disorder strength U the conduc¬ 
tance in the ribbon was clearly enhanced in comparison 
to the cylinder [see Fig. 2(b)], the conductance of the rib¬ 
bon for m > 0 is even smaller than in the cylinder [see 
Fig. 2(c)]. This behavior can be attributed to the absence 
of edge states at the sample edges for positive topolog¬ 
ical mass m > 0. In the cylindrical geometries we find 
that the delocalization transition is less pronounced for 
m > 0 than it was for m < 0 [compare blue dashed lines 
in Figs. 2(b) and (c)]. On the one hand the fact that the 
delocalization transition still exists for m > 0 supports 
our model of fiat bulk states undergoing a percolation 
transition. On the other hand, however, the more pro¬ 
nounced nature of the transition for m < 0 suggests that 
edge states propagating along the edges of the potential 
hills provide an additional link between localized states 
leading to a larger conductance. This picture, indeed, 
agrees very well with the analyses of Refs. 23 and 20, 
since the connecting local edge states in our model can 
directly be identified with the edge states that were pre¬ 
dicted to form in the non-trivial sub-gaps of the localized 
fiat bands. 

We would thus be in a perfect position to complement 
the theory of Refs. 20 and 23 with the intuitive expla¬ 
nation that these sub-gap edge states exist locally and 
connect bulk states localized around hills of the poten¬ 
tial landscape to form a percolating network of internal 
bulk states that lead to the decay of the TAI phase. The 
missing piece to complete our argument is to show that 
the picture we derived for the case of correlated disorder 
holds also for the uncorrelated case considered in Refs. 20 
and 23. 

We check this point explicitly, by verifying that our 
model can explain the appearance of the TAI as well 
as the observed delocalization-localization transition of 
the bulk states for the case of uncorrelated disorder. 
Consider, in this context, that the TAI conductance 
plateau in the ribbon geometry between U ~ SOmeV 
and U ~ 280 meV [see red curve in Fig. 2(a)] is destroyed 
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by the onset of the bulk delocalization at t/ 280 meV 
in the cylinder [see blue dashed curve Fig. 2(a)] which 
happens for much larger U than in the correlated case. 
Still, when the delocalization transition is in full effect (at 
U = 370 meV) the corresponding scattering states show a 
similar weighted distribution P{Ef — V) in the now spa¬ 
tially uncorrelated potential as already observed in the 
correlated case [see Fig. 5(d)]. Again the peak of this 
distribution fits nicely to the band structure of the clean 
limit, indicating that our picture of local edge states per¬ 
colating around internal edges of strong disorder holds 
also for the uncorrelated case. Last but not least, we 
mention that such a percolating state corresponds ex¬ 
actly to the ’’single bulk state” that is held responsible 
for the delocalization in Ref. 11. 


DISCUSSION 

Our results suggest that the emergence as well as the 
decay of the TAI phase depend strongly on the energy off¬ 
set and on the flatness of bulk valence bands in the clean 
limit. These flat bands feature an enhanced contribution 
to the density of states and occur in the center as well 
as at the boundaries of the BZ. Yet, the underlying BHZ 
model is only valid for small close to the F-point and 
thus does not yield a good approximation for the valence 
bands at the BZ boundaries of a real HgTe/CdTe quan¬ 
tum well.^ Correspondingly, we find that when chang¬ 
ing the grid spacing a in our discretized lattice from the 
value conventionally used in the literature (a = 5 nm) 
to different values, the position of the BZ boundaries 
= ±7r/a and the energy offset of these states at the 
BZ boundaries also change significantly. We also verified 
that the flatness of the bands at the BZ boundary is a di¬ 
rect consequence of the discreteness of the underlying lat¬ 
tice used for the numerical solution of the transport prob¬ 
lem (see also Refs. 6, 7, and 10 where discretized models 
were first employed to describe the TAI). As a result, the 
localization-delocalization transition and possibly even 
the TAI phase itself associated with these states at the 
BZ boundary will not occur in real HgTe/CdTe quantum 
wells as the strong-disorder limit in these devices will be 
different from the predictions of the discretized model. 
Quite remarkably, however, realizations of topological 
insulators have recently also been considered based on 
photonic systems.These so-called Floquet topologi¬ 
cal insulators are based on a discretized lattice of sites, 
just like in the numerical model used above to approx¬ 
imate the physics in HgTe/CdTe quantum wells. The 
strong-disorder physics, which we have discussed here, 
may thus well be realized in experiments based on effec¬ 
tive model systems in optics^^’^^ as well as in acoustics^^ 
or in other fields where wave scattering parameters can 
be tuned appropriately. 


CONCLUSION 

In this work we uncover the existence of percolating 
states in two-dimensional topological insulators. In par¬ 
ticular, we show how these states affect the phase bound¬ 
aries of the topological Anderson insulator, which is a 
topologically non-trivial phase caused by disorder. While 
the reason for the emergence of this phase has already 
been understood, its breakdown could so far only be 
vaguely connected to a delocalization of bulk states. Here 
we show that in a spatially correlated potential this de¬ 
localization is caused primarily by bulk states, that are 
localized when circumnavigating the hills of the disorder 
potential, but that become connected with each other 
when passing a percolation threshold. These connections 
and thus also the delocalization transition are consoli¬ 
dated by local edge states that can internally form in the 
disordered sample. By showing how the localized bulk 
states derive from flat bands in the valence band struc¬ 
ture of the clean sample without disorder, we clarify that 
the same physics is at work also in the well-studied case 
of an uncorrelated disorder potential. 


5 


uncorrelated 





correlated 



Figure 1. (Color online) Sketch of an uncorrelated disorder 
potential in a system of width and length W — L — 1000 nm 
and a correlated potential with correlation length ^ = 35 nm. 
The color code ranges from blue (strong negative potential) 
to red (strong positive potential). 



Figure 2. (Color online) Average conductance {G) as a 
function of disorder strength U through different HgTe/CdTe 
quantum wells. All systems are of width and length W = 
L = 1000 nm and the Fermi energy is set to Ep — 16meV. 
The ribbon geometries are shown by red solid curves while the 
cylinder geometries are shown by the blue dashed lines, (a) 
Average conductance through a system of topological mass 
m — — lOmeV (topological insulator) in an uncorrelated po¬ 
tential. The disorder average is taken over 1000 configura¬ 
tions. The ribbon geometry features the TAI conductance 
plateau of {G) — le^ jh for SOmeV ^ U < 280 meV. In 
this region the conductance through the cylinder is almost 
entirely suppressed followed by a delocalization transition of 
the bulk states starting at 280 meV. (b) Average conduc¬ 
tance through systems with the same parameters as in (a) but 
in a correlated potential with correlation length ^ = 35 nm. 
The cylinder geometry (blue dashed curve) clearly shows the 
bulk delocalization transition, which occurs here for weaker 
disorder strengths U than in the uncorrelated case, (c) Av¬ 
erage conductance through systems with positive topological 
mass m = lOmeV (ordinary insulator) in a correlated poten¬ 
tial of correlation length ^ = 35 nm. The disorder average is 
taken over 200 configurations. The cylinder geometry shows 
a less pronounced bulk delocalization transition than in (b), 
followed here also by the results for the ribbon due to the 
absence of edge states. 
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Figure 4. (Color online) Probability density distribu¬ 
tion P{V) for disorder values V{xi,yi) at grid points {xi,yi) 
weighted by the absolute value \'il;{xi,yi)\‘^ of the scattering 
wave function at Ep = IGmeV with injection in the first lead 
mode in a system of width and length W = L = 1000 nm and 
topological mass m = — lOmeV. A disorder average is taken 
over 1000 configurations of random potentials with disorder 
strength U = 220 meV and correlation length ^ = 35 nm. 
The distribution shows an enhancement of the wave func¬ 
tions at disorder values V situated between Vmin ~ 25meV 
and Fmax ~ 75meV. These values correspond nicely to the 
distance of Fermi energy Ep to the flat valence bulk band 
states at the BZ edges. 


Figure 3. (Color online) Comparison of the scattering wave 
function \'il;{xi,yi)\‘^ (shown on a logarithmic scale in the left 
column) and the potential itself (right column) at the delocal¬ 
ization transition {U = 220 meV) for systems with W — L — 
1000 nm, Ep = IGmeV and m = —lOmeV. We only colored 
the potential values between 26meV < V(xi,yi) < 78meV for 
which the Fermi energy is effectively shifted into the valence 
band of the clean band structure (the remaining potential val¬ 
ues are left in white). Two different disorder configurations for 
cylinder and ribbon geometry are considered in the top and 
bottom row, respectively. The similarity of the wave functions 
and these truncated potentials illustrates that the delocalizing 
bulk states are percolating around the hills of the potential 
landscape and that these percolating states have their origin 
in the bulk band structure of the clean sample. 
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Figure 5. (Color online) (a) Band structure E{kx) of a 
system of width W = 1000 nm and topological mass m = 
— lOmeV obtained by directly solving the Bloch-eigenvalue 
problem numerically for the Hamiltonian in Eq. (1). (b) Den¬ 
sity of states p{E) in a closed and clean system of width and 
length W — L — 1000 nm, topological mass m = —lOmeV 
and cylindrical geometry, (c) Weighted probability density 
distribution P(Ef — V) for systems with m = —lOmeV, 
W = L = 1000 nm in a correlated potential of correla¬ 
tion length ^ = 35 nm. The wave functions are calculated 
in the region of the bulk delocalization at disorder strength 
U = 220 meV and for 1000 random disorder configurations. 
The red solid and blue dashed lines show P{Ef — V) at 
Ef = IGmeV in the ribbon and the cylinder geometry. The 
green dash-dotted and black dotted lines represent P{Ef — V) 
for ribbon and cylinder at Ef = 30meV. The similarity of all 
four curves shows that the distribution is a very fundamen¬ 
tal system property. The shape of the distributions closely 
resembles the density of states at the BZ boundary shown in 
(b). This indicates the bulk-like nature of the states responsi¬ 
ble for the delocalization transition, (d) Weighted probability 
density distribution P(Ef — V) with the same system param¬ 
eters as in (c) but in an uncorrelated disorder potential. The 
pictures are again taken in the region of the bulk delocaliza¬ 
tion at disorder strength U = 370 meV and for 1000 disorder 
configurations. 
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